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The free energy and magnetization for the general SU(N) one impurity Kondo 
model in the magnetic field, h, are calculated by extending the previous 1/N expan- 
sion technique: the saddle point is determined self-consistently to the 1/N order. The 
obtained universal field dependent magnetization M(h/Tj() by this simple method 
is shown analytically to be asymptotically exact at both h <C Tr and h 3> Tk lim- 
its. For general "/-electron" fillings, except half filling, the M(H/Tk) curves cross 
continuously from weak to strong coupling limit, but overestimate the curvature in 
the crossover region for moderate N . The magnetic Wilson crossover numbers are 
calculated for amusement. Our results explicitly verify that the 1/N parameter is 
non-singular under the adiabatic continuation. 
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I. INTRODUCTION 



The flowing of an effective interaction from weak coupling at high energy to strong 
coupling at low energy is an important and frequently encountered phenomenon in various 
physical systems. A well known condensed matter example is the Kondo effect [?]. Usually 
it is only possible to construct perturbative solutions in the weak and strong coupling limits. 
That the Kondo problem admits an exact solution provides a useful testbed for new ideas and 
methods. Among various methods applied to the problem, the numerical renormalization 
group(NRG) [?], Bethe Ansatz [?], and probably Non-Crossing Approximation [?], nicely and 
accurately produce the crossover. Unfortunately, these methods either are very complicated 
or heavily rely on numerical calculations. A simple and elementary method describing the 
crossover is desirable and may give us new insight. 

Recently, motivated by the NRG results on the two impurity Kondo problem [?,?] which 
claim that there is a line of Fermi liquid fixed points continuously modified by the RKKY 
interaction between the two impurity spins, we have developed an "Eliashberg equation" 
approach to build the magnetic correlation between the two impurity spins nonperturbatively 
into the ground state [?]. Naturally, we want to test our method for the one impurity 
Kondo problem. In this simple case, our approach amounts to the self-consistent one-loop 
approximation. For the general SU (N) impurity spin model [?] with the orbital degeneracy 
N, we expand the free energy in 1/N and determine the saddle point self-consistently using 
the free energy including one-loop (1/iV) fluctuation contributions. We shall see that 1/N 
is a non-singular parameter under the adiabatic continuation [?], at least outside a narrow 
crossover region. The effect of high order terms is to smooth out the crossover. Technically, 
1/N fluctuations always involve cut-off dependent contributions. In order to obtain the 
universal free energy and magnetization, all the cut-off dependent terms have to be absorbed 
into the Kondo temperature Tk- In the following, we first sketch the procedure then give 
the details in the next two sections so that whoever not interested in details can skip from 
the end of introduction directly to the results. 
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The Kondo problem describes an impurity spin antiferromagnetically coupled with 
strength J to a wide conduction band with density of states p(e). The Hamiltonian for 
the general SU(N) model [?] in the magnetic field is 

H = £( CjE + oK)c\ c c- ka - L £ (clfMlc^,) + h £ a flf 9 . (1) 

k,a k,k',a,a' a=-S 

The impurity spin is represented by N = 2S + 1 localized degenerate levels partially filled 
with "/-electrons". Their creation and annihilation operators are subject to the constraint 

n f = J2f^ = %N. (2) 

a 

We have set the gyromagnetic ratio and Bohr magneton equal to one so that the magnetic 
field strength h has the energy scale. For Ce, the lower spin-orbit splitted multiplet usually 
has N = 6. The coefficient q is treated as a constant of order one [?] in the expansion and 
will be given any value at the end of calculation. We shall present results for q = 1/2 and 

go = l/N. 

There are two physical parameters in the Kondo problem, the bandwidth D and the 
dimensionless coupling constant g = Jp(0). In the scaling regime, h <C D and T K <C D, 
physical quantities depend on D and g only through the Kondo temperature T K = T K (D, g). 
If the initial bare g 1, we can find T K in the D/T K — > oo limit. This is equivalent to the 
ultraviolet renormalization. The renormalizability of the Kondo problem was stated long 
time ago [?,?] and can be proved without difficulty. After absorbing the bare parameters 
into T K , physical quantities such as the magnetization must be a one- variable function: 
M = M(h/Tjc), since M is dimensionless. Usually, there could be many different scaling 
functions M(x) with x = h/T K , depending on the band structures p(e) (cut-off schemes). 
However, M(x) for the Kondo problem is universal, because changing band structure only 
adds in irrelevant perturbations which quickly die out under scaling if initial g [?] . The 
only possible exception is particle-hole symmetry breaking perturbation which is marginal 
and may lead to a modified M(x). Thus, the obtained scaling solution for the magnetization 
in our calculation is directly comparable with any previous result up to a proportionality 
constant between different definitions of the Kondo temperature [?]. 
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It has been known from the phenomenology of dilute alloys [?] that the nature of 
the strong coupling fixed point of the Kondo problem is a local resonant level. The two 
parameters of the resonant level, its position e/ and width A, are precisely the saddle point 
parameters in the 1/N expansion [?]. Including 1/N fluctuations, the free energy in the 
magnetic field can be written as 



where the mean field and 1/N contributions, Fmf and F\m, have no explicit dependence 
on N. The two parameters e/ and A are determined by the stationary condition of the free 
energy. To find the Kondo temperature Tr-, we separate out from the free energy all terms 
depending on the bare parameters g and D, 



The regularized free energy, F reg , depends on g and D only through T K . With a proper defi- 
nition of Tk, F becomes a constant depending only on g and D, representing the correction 
to the ground state energy. The thermodynamics is contained in F reg from which we obtain 
the field dependent magnetization. 

The paper is organized as following. In the next section, we briefly recapture the large- N 
approach in the magnetic field to define our notations. The renormalization procedure is 
described in the third section. In the fourth section, we present the field dependent magne- 
tization from h <C Tk to h ^> T% for several values of N. The magnetic Wilson crossover 
numbers are calculated approximately. The proof that the magnetization calculated from 
F ieg has the correct h ^> Tk asymptotics and the integral expressions of some functions 
appearing in the regularization are included in the appendices for completeness. To alle- 
viate cross reference, we list the frequently occurring symbols together with their defining 
equation numbers in Table |. 



F(h, e f , A, g, D) = NF MF (h, e f , A, g, D) + F l/N (h, e f , A, g, D), 



(3) 



F(h, e f , A, g, D) = F(h, e f , A, g, D) + F Teg (h, e f , A, T K ). 



(4) 
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II. LARGE- N FORMALISM 



Following previous treatments [?,?], we introduce a Lagrange multiplier A to enforce the 
constraint @. By using the fact that the constraint commutes with the Hamiltonian, we 
write the partition function in the magnetic field h as 

(3dX 



^— Tr exp{-p[H + iX{n f - q N)]} 

Z7T 



V[c, c, f, /]exp 



- [ dr(C +H- iq N\) 
Jo 



(5) 
(6) 



After performing Hubbard- Stratonovich transformation to factorize the Kondo interaction, 
we rewrite the partition function as 

(3d\ 



Z 



2tt 



£>[c, c, f,f,Q,Q]exp 



P dr (Co + £' + ^j- - iq N\ 



(7) 
(8) 



The above Lagrangian possesses a U(l) gauge invariance 

/ ff -£ = / ff e* Q^Q = Qe~*\ A ^ A' = A 



dr 



(9) 



The redundant gauge degrees of freedom can be eliminated by choosing to work in the radial 
gauge. Separating the complex field Q into an amplitude and a phase Q = r e~ 1 ^, the phase 
<p can be absorbed into new variables f' a and A': f a = f a e - *^, A' = A + d<p/dr . In terms of 
new variables r, A', f a and f' a , the partition function can be cast in the form, after dropping 
the primes, 



Z = V[c, c, /, /, A, r] Y[ r(r) exp 



-J dTl£"(T) + — -tq NX 



(10) 



£" = E + H + + E /J(flr + «A + / i a)/ (T + E r + /< 



lcr C ko 



It is possible to completely gauge away the phase because it does not contain dy- 
namics. Since the last Lagrangian is bilinear in the Grassman variables c^ a and f a , we can 
integrate them out to obtain an effective action, 



Z = Z Q J V[\r\ exp[-S eff (A,r) + 5(0) £ dr\nr(r)] 
S cS = -J2 Tr In [d T + iX + ha + rG (r)r] + N dr (^j - iq X 

where 5(0) = (1/(5) J2„ n 1 with v n = 2vm//3, and 



(12) 
(13) 



Go(r) = -E 



(14) 



Z is the partition function of the non-interacting Fermi sea. 

The integration over the two real variables A and r can be expanded around a saddle 
point 



iX = 6f + iX, r = r + r. 



(15) 



Retaining only quadratic terms in A and r in the expansion, the partition function becomes, 
after dropping the tilde sign, 
Z 



= e -Sc ff (^o) f Y[ dX (u n )dr(u n ) exp 



(16) 



S$ = ^ EW-"n),r(-i/ B )) 



A(z/ n ) 



r(^ n ) 



(17) 



p(0)r 2 T x (u n ) ip(0)r T Xr (is n ) 
ip(0)r r X r(vn) p(0)T r (v n ) 

The zero temperature expressions of the matrix elements T's appearing in S^g have been 
given by Read [?]. Their extension to include magnetic field is straightforward. Here we 
have pulled out explicitly some prefactors for later convenience. 



F\(v n ) 

F\ r (v n ) 



E 



n Vkl(kl + 2A) 

2 



In 



4, + (kl + a) 2 



4. + A2 



E 



1 



Tr(^) = ^Ei 111 



" (f 
4 + (k| + A) 2 ] + 2A ln \e% + (k| + A) 22 



tan- 1 ■ e ?° 



(Tf) 



4. + a 2 



(18) 
(19) 
(20) 



where we have defined the mean field Kondo temperature, 



TP = Dexp[--\, g = Jp(0), 



1 



(21) 
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and the convenient notations, 



e f° = e / + a K 



7rp(0)r 2 . 



(22) 



The contributions to the free energy (Q) are given by 
1 



^tan- 1 ^ 

7T V A 



A l [ 4r + A 21 

— m — — 77\\ — 
2vr [ ( T f )2 



A (\ 



1 



'l/AT 



— E ln[r A (z/ n )r r (z/ n ) + r 2 r (z/„)] + const. 



(23) 
(24) 



In the free energy F\/n, we note that the prefactors in the front of T's in fll7Q exactly cancel 
the contribution Y^v n m r o °f originating from the Jacobian of transforming to the radial 



gauge. 



III. RENORMALIZATION 

To calculate zero temperature quantities, we can simply replace the discrete Matsubara 
frequency sum by an integration 

I pOO 1 pOO Ay 

F l/N = -j a^auv + iU W (25) 

The upper integration limit is actually cut off by the conduction electron bandwidth D. One 
can see this from the approximation we made in deriving the mean field free energy and 
1/N fluctuation matrix element T's, 

1 r D de f D \ 

E " = / • = -«2p(0) tan" 1 — ~ -Z7rp(0) sgnw n 9{D -\u n \). 

£ lUJ n - J-D tUJ n - 6 \U) n J 

Obviously, Fim of fl2"4|) contains contributions linear in D which become divergent in the 
D — > oo limit. A little investigation shows that the sub-leading divergent terms of -Fi/at have 
the form of In In D. 

To separate out the cutoff dependent terms of Fi/ N which diverge as D — > oo, we consider 
the v —>■ oo asymptotic behavior of the integrand, 



r» = r A r r + v\ r = - 

- 1£ 



i 

V 2 



^(lnz/) + ^r 2 (lnz/) + 0(u- 2 ) 



(26) 



The two functions T\ and only depend on In v and have the following simple forms 



T^lnz/) =4 

ryinz/) =4 



, 2 V , V 2/1 

In 7^ - ^2 In — + 7T - - g - ?7i 

J-K J-K \* 

A(l - tt?7 2 ) In ^- - tt^A Q - 7^ - vre/ 



(27) 
(28) 



where we have introduced following two short hand notations, 

_ dF MF (e f ,A) _ 1 1 r tfln -i^l 



dF MF {e f ,A) 1 
V2 = ET1T = T7Z 2^ ln 



4. + A2 ' 



<9A 



A% 



(o) 



(29) 
(30) 



K 



They both are independent of frequency v. The 1/N fluctuation free energy is regularized 
as following, 



lnTi(ln^) + 



2r 2 (lnz/) 



+ [ D ^InT^nv) + [ 



z/Ti(lnz/ 

u(/} rfx r 2 (x 



6{v - u ) 



+ const. 



(31) 



hn(u /T K ) TT Tl(x) 

Since the first integral is convergent, we have extended the upper integration limit to infinity. 
Note that is not a parameter of the theory. Fu^ is independent of uq. We shall choose 
it for computational convenience. Actually, it provides a useful consistency check for the 
numerical calculation. The cut-off dependence is then separated out from the last two 



integrals of (|3T|). 



1 f 

— / dz/lnri(lnz/) = £>Ai(£>, 7/1, rj 2 ) - v Ai(v , 771,772) 

,7T J vq 

I ,ln(D/T K ) Y 2 (x) 

dx - ; ; = A 2 (D, 771, 772, e/, A) - A 2 (z/ , 771, 772, e f , A). 



(32) 

7T J]n(uo/T K ) Tx(x) 

The so-defined two functions Ai and A 2 are given in the appendix. 

To treat the cut-off dependent terms D Ai(D, rji, 772) and A 2 (D, 771, 772, 6f, A), we first 
obtain explicitly 



A 2 (D,r) 1 ,r} 2 ,e f ,A) = — lnm — - 77 2 Alula—. 



(34) 
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where we have neglected terms which vanish as D — > oo. Using the fact that rjx and 772 are 
the derivatives of the mean field free energy, we can show that Ai and the second term of 
(|34|) can be renormalized away from the saddle point equations if we let the saddle point 
parameters e/ and A acquire following 1/N corrections, 

?/=«/ + ^A.M.>)2*) (35) 

A = A 4""4 + ^ Al(,: '* fl (36) 
where rjl and rj^ are the values at the point of the saddle point solution, e/ = e*j- and 
A = A*. When we rewrite the mean field free energy in terms of the renormalized saddle 
point parameters e/ and A, we have to include the difference F MF (ef, A) — F MF (ey, A) into 
the cut-off dependent part of the free energy F introduced in (f|). Collecting this term, 
(|34]), Ai(D,r]i,ri2), and a term coming from replacing T% by T K in F M f, the total cut-off 
dependent part of the free energy is 



A A D 

F = -N— In + — In In — + D 

rri 1 rri 1 

7T T K ix T K 



(37) 



Note that the last term is a constant, to the order 0(i]i) ~ 0(7] 2 ). The first two terms 
cancel out if we define 

Tk = Tf (\n^-Y 1/N = D fln^V^exp (-1) . (38) 



TkJ V T K J \ g / 

In the spirit of order by order renormalization, we replace ef, A and appearing in F 1 / N 
by 6f, A and T K respectively. This gives us the regularized free energy as a function of h, 
6f, A and Tk only. Note that our expression for the Kondo temperature is consistent with 
the well known expression Tk = Dg l / N exp(— 1/g) up to 0(1/ N). 

Actually, one can simply expand Ai(r/i, 772) in 1/N by using the fact r\\ ~ rj 2 ~ 0(1/N), 
a consequence of the saddle point equations. We immediately see that the only 0(1) contri- 
bution of Ai (771, 772) to the free energy is a constant. This constant is the correction to the 
ground state energy and has no effect on the physical quantities. Higher order terms in the 
expansion of Ai (771, 772) can be neglected in the order by order renormalization. The second 



term of ([34]) is also dropped since it is of order 0(1 /N). After we renormalize away the 
first term of (0) by defining the 1/N corrected Kondo temperature Tk via (|38| ) and replace 
the mean field Kondo temperature Tjp in T r by Tk, the resulting regularized free energy is 
then only a function of e/, A, h and Tk- All these are due to the fact that the free energy 
is stationary with respect to e/ and A. A 0(1/ N) shift of these parameters does not induce 
any change in the free energy to the order 0(N) + 0(1). 

After completing the renormalization, the universal free energy is, from 

A" A I ~ / v/ r ?„ + A 2 



and ©-(H, 



reg 



E 



^tan- 1 ^ 

7T V A 



7T 



N „ 



fo_ 

T 



K 



+N ^- - g j e f + Fyff 

F i e /N = ~ v Al(f , ?7l,772) - &2(u ,Vl,V2,£f,A) 

2r 2 (ln v) 



+ 



Z7T 



lnT^lnz/) + 



(1/ - i/ ) 

^(°) ic 



(39) 



(40) 



The parameters 7/1 and r] 2 only depend on e/, A. Inside r] 2 and T r , T£ J is replaced by Tr-. 

The saddle point parameters, e/ and A, are determined by solving the following two 
saddle point equations, 



^F re6 (M/,A,T x ) = i 



1 



an 



1 / 

A 



1 9 







1 (9 1 



7» 



A 2 



K 



d 

_ preg 

5A Viv 



0. 



(41) 
(42) 



Substituting the solution e/ = e*j(h/TK) and A = A*(/i/7V) back into F reg , we obtain the 
scaling form of the free energy depending only on H/Tk, up to an additive constant. The 
magnetization is 







M(h/T K ) = --^rF Ieg (h, e*f, A*, T K ) = ~£ a tan" 1 (-£ ) - ^ 



9/i 



7T 



d 



^reg 



a; ^ 



(43) 



The one-dimensional integration in the regularized 1/N free energy and its derivatives, as 
well as solving the two coupled equations ([41]) and (f42|) , are carried out numerically. 

We emphasize that the obtained magnetization is not a 1/N perturbative result if we 
solve the equations fl41|) and (|42| ) self-consistently, i.e. not by expanding e*j and A* in 1/N. 
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The fact that we only carried out perturbative ultraviolet renormalization only implies that 
the Kondo temperature defined by ( |38| ) is perturbatively accurate to the 1/N order. In other 
words, our result for F reg or M(H/Tk) is perturbative at high energy but not necessarily 
perturbative at low energy, depending on how we solve the saddle point equations. As we 
can see, the same renormalization procedure can be carried out for every physical quantity 
and their calculation is a straightforward exercise. 

IV. RESULTS 

The solution of the saddle point equations, e*f(h/T K ) and A*(h/T K ), for q = 1/6, iV = 6 
is shown in Fig. [I] as an example. Generally for go 7^ 1/2, there are more than one solution 
in the weak coupling regime for a given value of h/Tx- Certainly, the criterion is to choose 
one with the lowest energy. However, since we know the asymptotics at both weak and 
strong coupling limits, we can follow the solutions continuously by varying the magnetic 
field slightly each time. For q = 1/6 and N = 6 as an example, there are solutions other 
than that shown in Fig. |l]for h/Tx > 0.52 and give magnetizations much closer to Hewson 
and Rasul's exact results [?] in near crossover region compared with the results shown in 
Fig. H But, if we follow these solutions to high magnetic field, they do not have the correct 
asymptotics. 

The field dependent magnetizations M(h/Tjc) for g = 1/2 and various values of N are 
shown in Fig. |2|. Note that each curve has a window in the crossover region where no solution 
is found by the present method. This happens only for go = 1/2- The reason is following. 
We try to describe the strong coupling fixed point by a resonant level. The particle-hole 
symmetry presented in the go = 1/2 case ties the position of the resonant level at the Fermi 
surface, e*j = 0, in the strong coupling regime. Certainly, the nature of the weak coupling 
is no longer a resonant level, thus e£ 7^ 0. A discontinuity must occur at some value of 
with increasing magnetic field h, preventing continuous crossover from one side to the other. 
Nevertheless, the window quickly narrows with increasing N. For N = 8, the solid line of 
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Fig. 0, the window narrows to 0.45 < h/Tx < 0.55. The indication is that probably we 
need infinite order of terms in 1/N to close the window and to obtain completely smooth 
crossover. The more terms we put in, the better the quality is in the crossover region. Similar 
features can also be seen for general values of go- in Fig. |3|, we show the magnetizations 
for q = 1/N, the "realistic" situation. Also shown are Hewson and Rasul's Bethe Ansatz 
results [?,?] for N = 6, 8. Although the lines can cross continuously from one side to the 
other, they obviously overestimate the curvature in the crossover region. With increasing 
N, the curvature is reduced. 

For amusement, we calculate the magnetic Wilson crossover numbers for the Coqblin- 
Schrieffer model [?], q = 1/N, although the calculation can be done for other values of q . 
The ambiguity in relating Tk from different cutoff schemes can be eliminated by imposing 
the condition of a vanishing hi~ 2 {h/Tx) term in the H/Tk 3> 1 expansion of M(H/Tk). The 
weak coupling scaling form for the magnetization in terms of Tk is well known [?], 
M 1 In In ^- In 2 

= 1 z t e tt + 5-E- + -", h/T K ^l. (44) 

M 21n#- 2iVln 2 ^ N\n 2 ^- 1 K J 

±k Ik Ik 

The last term of (|4"1) can be removed by changing to a new energy scale 

21n2\ 



T h = 2-V N T K ~T K / (l 



(45) 



N J 

Although we only explicitly prove the first log term of fl44|) in the appendix, we expect 
that our result ( |l~3|) will precisely produce all three log terms of (|44]), since all 1/N order 
contributions to the free energy are included in the present approach. Another direct way to 
see this is following. Given the second term of (f4"4"D, the last two terms of (44) are determined 
by the second term of the weak coupling beta function [?] , 

Our expression for the Kondo temperature fl5!f ) gives exactly the same beta function. The 
correct asymptotic form ( |P| ) allows unambiguous determination of the energy scale Th in 
the present approach. In terms of the unique energy scale Th, the coefficient a' in the strong 
coupling asymptotic form of the magnetization 
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is just the magnetic Wilson crossover number. From (f|7|) and (j45[) , we see a' = a/(l + 
21n2/iV). The slope a will be determined directly from M(h/Tx) curve. We list the results 
for the general SU(N) cases in Table ||. 

In summary, we calculated the universal field dependent magnetization for the general 
SU(N) one impurity Kondo model for various values of N and "/-electron" fillings. At both 
small and high field limits, our results become asymptotically exact, as shown analytically 
in the appendix. For other than half filling of the " /-electrons" , the magnetization curves 
cross continuously from one side to the other. In the crossover region, the bigger is the N, 
the smoother and the more accurate is the magnetization. In contrast to a continuous phase 
transition, the crossover involves no divergence. The other facet of the story is that one then 
does need high order terms to smooth out the crossover for a given N. 

It is my pleasure to thank Ian Affleck, Natan Andrei, Piers Coleman, Eric Sorenson, 
Eugene Wong for many helpful discussions. Illuminating discussion with Affleck helped me 
to understand many crucial points. I am also grateful to Alex Hewson. He kindly provided 
me the Bethe Ansatz results. This work was supported by NSERC of Canada. 

APPENDIX A: INTEGRALS Ai AND A 2 

For simplicity, we set T# = 1 in this section. From the definition, A 2 is an integral of 
the type, 



A 2 (D, T]!, t] 2 , e f , A) - A 2 0o, Vx, V2, e/, A) = 



D dv 



w In v + v 



(Al) 



In 2 is + a\nv + b 



where a, 6, w and v are all independent of frequency and are given by 



a(e f , A) = -7i7] 2 



(A2) 




(A3) 



w(e f ,A) = A(1-tt77 2 ) 



(A4) 
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*>(e/, A) = -vrAr/s - vrr/ 2 ^ - 7r e/ f- - q - r)t ) . 



(AS) 



By carrying out integration, we find 



A2O0, V11V2, e f , A) = ^- In (in 2 z/ + a In u + b) 



2v-aw J ^gln( ;!::!gg ), a 2 -46>0 



2tt 



a 2 - 46 < 



(A6) 



— In (in 2 z/ + a In v + 6 



(A7) 



VifoZ^ """" V21n^o+a/ ' 

From the definition of Ai, it is an integral of the type 

' D dv 

where we choose vq big enough so that the argument of the log function is always positive. 
We can see that Ai(uq, 771, 7/2) is analytic in a and b for small values of a and b. In some cases, 
Ai can be expressed in terms of the standard integral of exponential functions such as Ei(x). 
In the present problem, the parameters a and b never get very big. A series expansion is 
sufficient for the practical purpose. The expression we used in the present calculation is, 

'P 1 (\n- 1 u ) Ei(\nv o y 



7tAi(i/ , 771,772) 



In i/ 



+ E(-i)" +1 ^^K + ^) + 2 



eV /a2/4-b + e -y/a?/4-bJ g -a/2 



n=l 



In In vq 



In i/ 



where P n are polynomials of In 1 1/ , 



(A8) 



p n (x) = 1 + nx + n(n + l)x 2 H h n(n + 1) ■ ■ • (m - l)x m ~ n , 



(A9) 



and a%, a 2 are related to a, b through 



a± + «2 = a, 



ai«2 = 6- 



(A10) 



Ei(x) is the standard integral of exponential function, defined by 



Ei(x) 



—00 t 



-dt. 



(All) 



Note that a" + are expressed as polynomials of a and b. In the expansion ( |A8| ) , m is 
the order of expansion. The neglected terms are of the order [Max(|ai|, |a2|)/hiz/o] m+1 / m - 
Typical values used in our calculation are m ~ 10 — 15 and In uq ~ 5 — 8. 
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APPENDIX B: HIGH FIELD ASYMPTOTICS OF THE MAGNETIZATION 



The small field asymptotic behavior of ( fT7| ) is the well known result of the present 
approach [?]. Here, we prove the high field asymptotics for go = 1/2 and go = 1/-ZV. The 
proof for other values of go goes parallel. We shall set Tk = 1 and omit the star sign in the 
notation of saddle point solution e*f(h) and A*(h). 

Let's first consider g = 1/2 and even N. In the high magnetic field, the "/-electron" 
level is split into N levels. Each of them is distant from the others. For g = 1/2, the 
"/-electrons" occupy the lowest N/2 levels: a = —S, —S + 1, ■ ■ • , —1/2. The o = — 1/2 
level will lie close to the Fermi level. Spin exchange will result in a small resonant width. 
Thus, we write the solution in the form, 

h 



-6€i 



6 Jl A 

h ' h 



0. 



as h 



oo. 



(Bl) 



We recall that S is the spin and N = 2S + 1. Since we are looking for In 1 h asymptotic 
terms, we neglect all terms which die as h~ x or faster. Thus, 



€f + ah 



a 



h. 



-5e 



a 



(B2) 



With this approximation, the magnetization is simplified to 



M = Mi 



o 



1 

2^ 



tan 



Set 



D dv 1 



o 2-nYiv) 



v ( \ dT r^r ( \ dT ^or ( ^ 



dh 



dh 



dh 



(B3) 



where M = J2cr>o a i * s ^ ne saturation value of the magnetization. To shorten the notation, 
we use the unregularized 1/N fluctuation energy (^) to carry out the proof. Since the 
values for 5ef and A are given by the saddle point equations fl4lD and (fl||), we have to make 
use of them. With the simplification (|B2j ), The equation fl4lf) is similarly reduced to 

dT Xr ' 



1 i / A \ r D dv 1 

•-tan — + / r~, 

7T \6e f ) Jo 27rr(z/) 



,dT r dT x 
— h r r (z/) 



2T Xr (u 



de f ~'^'de f 

The matrix element T's involve the spin component summation J2 a 



de, 







(B4) 



1 



nV Tx 



V) 



N 



Ar 1 



V) 
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Each spin component of the T's can be read off from (|l8l)- (|20|) . The difference between 



the derivatives of the 1/N free energy appearing in ( p3|) and (|B4|) is that d/dh in ( P33|) 
will bring down an additional factor o with respect to d/dej. Dividing ( |B4|) by two and 
subtracting it from (p3|), we find 



+ TJv) 



dY 



+ 2YM 



dY 



Xr 



• (B5) 



9e/ 9e/ de/ 

Note that the a = —1/2 component vanishes in the above a summation so we can replace 
e /o- by (o + l/2)h. Carrying out the derivatives, we find 

h(o + \) 2 \ 2 



M = M, 



o 



1 y, r 1 ' <h- 



txN 



Jo I» (a + ±) 2 /i 2 + + A) 2 \N^ s ln 



e% + + Af 



A 2 



v 



4(i/ + A) 



Nh(a + 



2/ A* 



tan 



-i 



6 /m 
i/+ A 



- tan" 1 f % 
A 



• (B6) 



iV(z/ + A) 
By noting, from the equation (|42|) , 

1 E ln(4 + A 2 )~(9(l/iV), 

fj,=—S 

we can expand the expression inside curly bracket of (|B6|) in 1/N. We shall also expand 



r(v) 



^ E in(4 + (- + Ay 

v it=-S 



tan 



6 /M 

+ A 



— tan 



v+ A 



+ 0(1/JV). 



(B7) 



By changing the dummy variable, v — hx, we can make following expansion 
^ E ^ [4 + (" + A )1 = 21n/i + ^ E \n{(S + rf + x 2 

[M=—S fl=—S 

= 2 In /i [1 + (In a;/ ln/i)] , 

where we dropped terms of order A/h as usual. That it is possible to make lnx/ln/i 
expansion in the last expression is due to the convergence of the integration in (|B6p. We 
also expand Y(u), given by ([B7|), in In -1 h and keep the leading term. The upper integration 
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limit in (|B6| ) can be extended to infinity. The final result for the magnetization is, after 
some manipulations, 

*L = 1—±Y, r~ {a + f -i-[l + Q(lnx/lnfe)] 

A#_ \r»yf_ ^ L „ / i\2 „ l„ J, L V / )\ 



Mo NMo^J, ir (^ + iy + x 2\uh 

1 1 ^ , 1. 



E 



(7 



In 2iVM f^, 1 2 1 

= 1 -FL + ° (ln " 2ft) ' (B8) 

For go — l/-^> strictly speaking, 1/N is no longer the loop expansion parameter. Never- 
theless, if we repeat the above steps, we find 

Note that the leading log correction is independent of N for q = 1/N. It is easy to see 
this from the perturbation in g. This term comes from the linear term, g/2, in the g < 1 
perturbation. The diagram for this term involves one conduction electron loop and one 
"/-electron" loop which together contribute a factor N 2 . The interaction vertex brings in a 
factor 1/N. After normalization, i.e. dividing by S ~ N, it is independent of N. 
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TABLES 



TABLE I. Definition of symbols and notations 



oymbol 


Definition (Eq. No.) 


symbol 


Definition (Eq. No.) 


D 


Bandwidth 


L 


© 


P(e) 


Density of states 


i i 

Li 


/lr) r-|\ 

(EO) 


h 


® 


L 2 






(0) 


'/i 


/lnnl\ 


e/, r 


(y) 


'/2 


/1orJ\ 

(|30|) 


La 




A-i 


Q32|), (|A8D 


i 

1 Ar 




A 

A 2 


/looh ifjlh 4X73"! 

(1221), (e^), (|A^) 


L r 


(H) 




© 


o T (0) 


(!) 


T K 








F 


(I), © 


Fmf 


©, (ED 


Freg 


(D, dD 


Fl/ N 


©, (H 


*1/N 





TABLE II. The calculated magnetic Wilson crossover numbers for the Coqblin-Schrieffer 



model, qo = 1/-ZV, defined as a' of (|47|). With Tr- defined by ([381), we read off the initial gra- 
dient, a in (|43|), the magnetization curve. Then the crossover number is a' = ct/(l + 2ln2/N). 

N Crossover number Bethe Ansatz 

2 0.25 0.342 (=1/^) 

4 0.65 

6 1.01 

8 1.36 

10 1.70 
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FIGURES 

FIG. 1. The solution of the saddle point equations in the magnetic field for qo = 1/6 and 
N = 6. Tk is defined by (|38|). e/ is the position of the resonant level and A is the width. 

FIG. 2. The universal magnetic field dependent magnetization for go = 1/2 and for N = 2 (short 
dashed line), N = 4 (long dashed line), N = 6 (dash-dotted line), N = 8 (solid line). All curves are 
parameter free. Note the improving quality for larger N . 

FIG. 3. The universal magnetic field dependent magnetization for the Coqblin-Schrieffer model, 
i.e. q = 1/N, and for N = 6(dashed line), N = 8(dash-dotted line), N = 10(solid line). All 
curves are parameter free. The points are Hewson and Rasul's Bethe Ansatz results: ./V = 6(filled 
triangles), N = 8(filled circles). The proportionality factor between Tk defined by ( |3"8| ) and the T\ 
appearing in Bethe Ansatz solution is determined for each N by matching the small field gradient 
of the magnetization. 
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